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ABSTRACT 


In the late 1970 *s, the Maryland National Capital Park and Planning Com- 
mission, the county level planning offices for Montgomery County, Maryland, 
was faced with the problems of managing the rapid growth of their 400 
square mile jurisdiction. The planning board decided digital data bases 
would be constructed capable of site planning and defining parameters nec- 
essary for energy, planning, hydrologic, economic and various resource for- 
casting models. Techniques for managing and collecting regional ancillary 
data (aerial photos, soils, etc.) and digital remote sensed data sets would 
be required to run models and produce timely graphics/statistics for 
present/proposed decision-making. Concurrently, the University of Maryland 
Civil Engineering Department (UOMCE) was developing a computer-based 
gridded information system (GIS) allowing engineers to create, access, 
integrate, and maintain a multi-parameter geographical data base for real- 
time hydrologic modeling. 

By joining forces, the UOMCE and MNCPPC developed an interactive/batch GIS 
(array of cells georeferenced to USGS quad sheets) and interfacing applica- 
tion programs (e.g., hydrologic models). This system allows non-programmer 
users to request any data set(s) stored in the MNCPPC data base by inputing 
any random polygon's (watershed, political zone) boundary points. The data 
base information contained within this polygon can be used to produce maps, 
statistics, and define model parameters for the area. Present/proposed 
conditions for the area may be compared by inputing future usage (land 
cover, soils, slope, etc.). 

This system, known as the Hydrologic Analysis Program (HAP), is currently 
operational on the MNCPPC 's HP 3000 mini-computer and the UOMCE UNIVAC-1180 
main-frame computer. HAP has been especially effective in the real-time 
analysis of proposed land cover changes on runoff hydrographs and 
graphics/statistics resource inventories of random study area/watersheds. 
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INTRODUCTION 


In August of 1977, the Montgomery County Offices of The Maryland 
National Capital Park and Planning Commission incorporated an adaptation 
of a computer-based geographical information system known as MSDAMP. 

MSDAMP, an acronym for Multi-Spatial Data Analysis Mapping Program, was 
developed at the Iowa ,State University Land Analysis Laboratory in November 
1972. This effort allowed cultural and physical data from existing maps 
and aerial photographs to be converted to a digital format and stored on 
the County's computer. The dominant land use or other desired parameter 
for each of these cells is encoded and entered into the computer. In batch 
mode, MSDAMP can be used to produce line printer gray-scale maps showing 
the distribution of various land uses, geologic features, slopes, soils, 
etc. 


The data structure of MSDAMP is an array of five-second cells covering 
the entire 400 square mile Montgomery County. Each cell, encompassing an 
area of 4.58 acres, measures 397.75 feet east/west and 505.90 feet north/ 
south. These dimensions were selected to give the five second increment at 
a latitutde of 39°00*00" and allow distortion-free symbolic maps to be pro- 
duced with ten column by eight line high speed line printers. At this cell 
resolution, each data plane (land cover, soils, etc.) in the Montgomery 
County data base would contain 55895 cells. MSDAMP requires the input of 
each individual cell georeferenced by latitude and longitude. In the sub- 
sequent years MSDAMP proved to be generally ineffective because of this 
cumbersome data collection method, inability to interface with application 
programs (e.g., hydrologic and planning models) and inefficient sequential 
data base searches. 

In 1978, the University of Maryland Civil Engineering Department 
(UOMCE) was conducting research in geographical information systems (CIS) 
for hydrologic analysis. The goal of the UC^CE CIS was to improve large 
area geoencoding techniques, provide a data structure capable of manipulat- 
ing random polygons of data within the data base, interface the data base 
with existing hydrologic models, and manage the data base consisting of 
both digital remote sensed data (Landsat, digital terrain) and ancillary 
hard copy information (aerial photography and soil maps). 

By joining forces, the MNCPPC and UCWCE created a system called the 
Hydrologic Analysis Program (HAP) which satisfied both organizations goals. 
The MNCPPC involvement ensured that the HAP software/hardware requirements 
would be developed for a non -programmer county level production environment. 

It was decided that the HAP data collection, data management, and 
model interface would be tested by expanding the MNCPPC data base to in- 
clude land cover, hydrologic soil groups and slope necessary for defining 
parameters in the Soil Conservation Service hydrologic model, TR-55. In a 
typical real-time situation, the planner types the coordinates of a water- 
shed boundary on the keyboard of an office terminal connected with a 
Hewlett-Packard 3000 mini -computer. The appropriate HAP compatible MSDAMP 
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data files are automatically accessed and the planner's terminal immediate- 
ly produces: 

1. a series of symbolic maps showing the distribution of land use, 
soil type, and slope; 

2. statistical tables showing the number of acres and the percentage 
of the watershed area devoted to each land use, soil type, or 
slope category; 

3. a table showing the runoff curve number, hydraulic length, and 
time of concentration needed to enter the SCS model. 

If the planner also enters a rainfall amount into the terminal, HAP 
will automatically list the volume of runoff and estimated peak discharge 
as computed with SCS-TR-55. 

A major function of HAP is to give planners the capability to assess 
the impact of land use or other changes being considered for the watershed. 
Thus, after the above information has been obtained for existing conditions, 
the planner can type the coordinates of those subareas being considered for 
change into the keyboard. He lists the type of new land use or other 
changes being considered for each of these subareas. The program then re- 
cycles in the manner outlined above with the new subarea information to 
give data for the watershed. under new conditions. 

HAP MODEL 


The guiding principle in the development of HAP was that planners must 
be able to obtain quantitative information in real time for any watershed/ 
political zone in the County. This pilot HAP model consists of three major 
components: 1) an enhanced CIS designed for collecting and manipulating 

the data base; 2) the interface between the data base and SCS-TR-55 hydro- 
logic model; and 3) graphics/statistic package. 

HAP Gridded Information System 


The HAP CIS was designed to simplify the user's maneuvering of data 
between the hierarchical levels of the data base. This task was accom- 
plished by providing the user with a hierarchical STUDY AREA-7J$° UNITED 
STATES GEOLOGICAL SURVEY (USGS) TOPO SHEET-GRID-CELL-ATRIBUTE model. 

Figure 1 shows the Montgomery County, Maryland study area subdivided into 
its nineteen USGS topo sheets. Each topo consists of 90x90 grid of 4.58 
acre cells containing various geophysical attributes. The cells are geo- 
referenced by latitude/longitude and grid row/col. 

The concept of storing geographical information in a computer can be 
illustrated by Figure 2. Conceptually, the spatial distribution of the 
geographical quantities are coded as an array of rectangular cells with the 
position of each cell identified by a row and column number. The dominant 
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geographical quantity within each cell is identified, generally, by a 
single valued alphanumeric character. This information is stored in the 
computer in some format similar to that listed in the lower part of Figure 
2. The strategy illustrated by Figure 2 can be extended to a system in 
which layers of geographic information are reduced to the same format, 
added to the computer and interfaced with each other to produce multi- 
variant parameters, such as those required for hydrologic models. Figure 3 
is a schematic illustrating the organization of a multi-parameter data base. 

Using the HAP model as a loading template, HAP will accept data by the 
gird-overlay featured in Figure 2, digitized polygons or previously existing 
gridded digital data. The electronic digitization of single-valued polygons 
is very efficient for homogeneous data attributes like hydrologic soil 
groups. The polygons are converted to grid cells within HAP. To minimize 
collection efforts, an entire topos grid is tagged with the dominant data 
class. Only the cells containing subsequent data classes are actually geo- 
encoded. 

The model also allows an engineer or planner to use the quad sheet as 
a ’’pointer” and, thereby, manipulate only the data within the quad sheets 
involved in an analysis, rather than the entire MSDAMP data base. For 
example, if a planner was interested in the analysis of the shaded watershed 
shown in Figure 1, the first step would be to type in the topo access num- 
bers 2542, 2819, 3650, and 3927. This would cause the computer to access 
the tape or other off line inexpensive storage and bring the 90 x 90, 5- 
second arrays of cells contained within the Gaithersburg, Sandy Spring, 
Rockville, and Kensington quadrangle sheets into temporary direct access 
storage. The role of the topo access will be explained in a subsequent 
section. 

Table I lists the land use categories accessed by HAP from the MSDAMP 
data base. Unless specified by the planner from the terminal, the symbols 
listed in Table I will be used by HAP in printing any maps requested in 
the output. The planner has an array of options including the assignment 
of ”blanks” to some of the categories, or all but one, in order to produce 
special purpose thematic maps. The Curve Numbers listed in Table I will be 
discussed in the following section. 

Table II lists the slope categories and their map symbols stored in the 
HAP data base. Ranges of slopes, rather than specific slope values, are 
used in order to allow the slope within a cell to be represented as a single 
digit and, thereby, minimize the storage requirements. 

HAP Hydrologic Model 

The power of the system is realized when inese arrays of stored vari- 
ables are interfaced directly through the computer with simulation models. 

In this approach, cells within a watershed are combined to define the input 
parameters needed for the hydrologic model. The model then outputs desired 
streamflow characteristics in a format appropriate to the user’s require- 
ments . 
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Used In Hydrologic Analysis Program 
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Table 11 

Slope Category Symbols Used In Hydrologic Analysis Program 
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The model used to {jenerate the volumes of runoff and peak discharges 
from the information available in the data base is a computerized version 
of parts of SCS-TR-55^. It is assumed that the HAP user has or will develop 
an understanding of SCS-TR-55. Thus, only the key equations needed to ex- 
plain the logic flow of HAP will be presented. 

The implementation of HAP starts with the user drawing a boundary 
around the watershed of interest on USGS minute quad sheets. Names of 
the quad sheets involved and a sufficient number of points to adequately 
define the boundary are entered into the terminal. The internal software 
of HAP then isolates the cells contained within the watershed boundary. 

The land use and hydrologic soil group are interfaced to compute a curve 
number for each cell in accordance with Table I. The cells are then 
summed and divided by the total to obtain an average Curve Number for the 
watershed. In a similar fashion, the average slope is obtained for the 
watershed. Table III lists the symbols that will be used to print HAP 
soil maps if desired. 

After the Curve Number, CN, has been derived, the relationship 


S = 


.1000 

^ CN 


) - 10 


( 1 ) 


is then used to compute the potential maximum storage, S, of the watershed. 
The result of this computation is then entered into 


Q = 


(P-.2S)^ 

P+.8S 


( 2 ) 


to obtain the volume of runoff, Q, from the rainfall, P. 


HAP obtains the area of the watershed. A, by counting the mamber of 
cells encompassed within the boundary and multiplying by 4.58. The area 
is then available to estimate the hydraulic length in accordance with 

= 209 (A) (3) 

The time of concentration is then estimated from 


T 

c 


1 

0.6 


• 8 >7 

H “(S+1) 

1900 Y 


(4) 


where Y is the average slope of the watershed. Finally, the time of con- 
centration is entered into a mathematical relationship defining the curve 
of Figure 4, to produce a dimensionless peak discharge which HAP converts 
to the peak discharge in cubic feet per second by multiplying the area in 
sq. miles by the volume of runoff obtained from Equation 2. 
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The use of Figure 4 imposes the requirement that the rainfall used 
with HAP be a 24 hour volume so it will be consistent with the intensity 
distribution of the SCS Type-II Storm Rainfall. The SCS Type-Il Storm is 
excellent for use in urban and urban fringe studies. The storm structure 
provides a period of several hours of light rainfall during which the soils 
are "wetted-up" to reduce infiltration rates and fill a portion of the 
depression storage. There is then a period of several hours of heavy rain- 
fall simulating thunderstorm intensities. Finally, there is a period of 
lightening rainfall during the last several hours. 

It should be recognized that there are a number of options within SCS- 
TR-55. One approach within TR-55 is to use detailed maps to measure in- 
cremental hydraulic flow lengths and then develop the time of concentration 
through velocity computations. HAP defines the hydraulic length and the 
time of concentration with the empirical Equations 3 and 4 which were 
developed by SCS from regression analyses of watersheds throughout the 
United States. A constraint on Equation 4 is that it be limited to water- 
sheds of less than 2000 acres. Table IV contains Montgomery County rainfall - 
frequency values required by TR-55. 

Table IV 


24 Hour Rainfall in Inches for Different Frequency Events 

for Montgomery County 


Frequency 


Rainfall Inches 


1 

2 
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10 

25 

50 

100 


2.6 

3.2 
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5.1 
5.6 

6.3 

7.2 
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USE OF THE HYDROLOGIC ANALYSTS PROGRAM 
FOR WATERSHED INVESTIGATIONS 


The use of HAP as a tool for watershed investigations will be illus- 
trated through the use of an example problem. The specific problem ad- 
dressed is the comparison of the peak discharge for a watershed under 
present conditions with the peak discharge anticipated under conditions of 
complete development as allowed by the Montgomery County approved land use 
master plan. The statement of the example problem is as follows: 


There is a need to provide stormwater storage behind and existing 
conduit on the Rockville quad. The design of the control struct- 
ure that will create the backwater requires estimates of the peak 
discharge for present and future land use conditions for a 100 
year, 7.2 inch, 24 hour storm. Determine the current and future 
peak discharges, volumes of runoff and percent changes created in 
these quantities by developing the watershed to its ultimate land 
use distribution. 


HAP uses four "modules” during its execution. The function of these 
modules are: 


BOUNDARY MODULE* The BOUNDARY module is called when the user wishes 
to enter the coordinates of the watershed boundary. The BOUNDARY 
module retrieves and stores the cells contained within the boundary 
of the study area. 

HYDROLOGIC MODEL* The HYDROLOGIC MODEL module merges the cells 
within the boundary and computes the slope for the study area. 

These parameters are then interfaced with the SCS-TR-55 and 
volumes of runoff and peak discharges estimated for user input 
rainfalls. The user is also prompted for map production and 
statistical summaries of the study area for present conditions 
(i.e., land cover, soils, and slope). 

UPDATE MODULE* The UPDATE module allows the input of proposed cell 
changes within the study area. The cell distribution in the water- 
shed is adjusted for the proposed conditions interfaced with 
TR-55 and revised runoffs compared to the present conditions. 

The user is prompted for production of maps and statistical 
summaries describing the study area under proposed conditions. 

STOP MODULE* The STOP module is called when the user has completed 
his study and wishes to exit from the HAP program. 
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Hie solution of the example problem will assume the user is accessing the 
program capabilities with a 132 character terminal. The program is interactive 
with the terminal outputing prompts and explanations to aid the user in inputing 
data or requesting HAP operations. Those terminal statements within boxes are 
input statements entered by the user. Those statements in upper case letters are 
outputs from the terminal. 

The following run stream describing the steps in the solution of the example 
problem contains explanations of key points. 


TYPICAL HYDROLOGIC ANALYSIS PROGRAM 
RUN STREAM 


STEP ONE GATHER NEEDED INFORMATION PRIOR TO ACCESSING HAP 

Before working with the terminal, it is necessary, in the interest 
of efficiency, to tabulate the information necessary to define the 
watershed boundary and move the appropriate portion of the data 
base into direct access storage. The first step is to assemble 
the topographic sheets needed to define the watershed boundary. 

The watershed boundary is sketched. The 90 x 90 transparent mylar 
grid is overlaid on the quadsheet. Enough nodes are picked to allow 
the watershed boundary to be approximated as a series of straight line 
segments. In this example, the entire watershed is within the bounds 
of the Rockville quadrangle sheet. In order to get the data of this 
quadrangle into direct access storage, the user will have to enter 
the disk access number obtained from Figure 1, in this case 2819. The 
user will also have to enter the latitude and longitude of the upper 
left hand comer of the quadrangle sheet. Finally, the cells to be 
changed to a proposed land cover will have to be tabulated for entry 
from the terminal. Thus, prior to accessing the terminal, the user 
would write down the information below. 

Utitude » 390730 
Longitude * 771500 
Disk Access No. * 2819 


193 


Tabulation of cell coordination of watershed boundary. Proceed clockwise. 
The boundary will automatically be closed between the first and last points 
listed. 


Watershed Boundary 


Boundary Point Bov Column 
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Tabulation of cells to be changed to reflect proposed watershed 
conditions. The cell location and proposed land cover change within 
each cell is tabulated using the transparent mylar grid and the land 
cover codes of Table I . If soil type or slope were to also be changed, 
Tables II and III would be used to obtain the proper symbols. 
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If more than one quadrangle is involved, this process will 
be repeated for each sheet. 
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STEP TWO SIGN ON HEWLETT PACKARD 3000 AS FOLLOWS 


(a) Turn on Terminal 

(b) Wait for green light 

(c) Hit return key 

(d) You will get a colon back 

(e) Type in (HELLO USER. MPLENV, WATER) 

(f) It will ask you for Password 

(g) Type in ENVIRON 

(h) T^e in NAZ 

(i) You will get # sign 

(j) To execute HAP you type in HAP 


THE TERMINAL INPUTS AND OUTPUTS FOR STEP 2 ARE: 


MlHtfiHttititifiHiHYDROLOGIC ANALYSIS PROGRARMtuiiM 
««£NTER: BOUNDASY.HYMOLOGIC NODEL.UPOATE.LOAO.STOPm 


Note: At this point the terminal will stop. The output immediately 

above lists the modules available. HAP is ready to accept the 
input of the watershed boundary. 


STEP THREE BOUNDARY MODULE 


A. Request the BOUNDARY module by. typing BOUNDARY into the terminal. 

B. Specify if you need input prompts. 

C. Enter name of study area. 

D. Enter number of topos containing study area. 

E. Enter latitude and longitude of northwest topo comer. 

F. Enter DAF topo address (see Figure 1) 

G. enter whether boundary points are P=POLYGON or C=CELL format. 

H. Enter whether boundary points are I=INTEGER or D=DIGITAL. 

I. In this example the boundary is defined from the northeast 
comer row and column coordinates of a 90 x 90 (4.58 acre) 
cell grid. If the boundary points were entered from a digit- 
izier the map scale would be requested. The program is set 
up to accept the node points from a digitizer in inches 
(eight sets of X, Y coordinates per card in F5.2 format). 
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TOE TERMINAL INPUTS AND OUTPUTS FOR STEP 3 ARE: 


<Mia«nrl 
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STEP FOUR HYDROLOGIC ANALYSIS 


The Hydrologic Analysis Program heading listing the modules is on the terminal 
output. The cells within the watershed boundary have been assembled. HAP 
is ready to do the TR-55 analysis. 

A. Request the HYDROLOGIC MODEL module. 

B. Input rainfall (inches). 

C. Specify if present condition land cover, soils, and slope 
maps are to be generated and to what unit they should be 
transmitted. 


TOE TERMINAL INPUTS AND OUTPUTS FOR STEP 4 ARE: 
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STEP FIVE 


ENTER PROPOSED CHANGES OF STUDY AREA AND COMPARE SCS TR-55 
RESPONSE FOR PRESENT/PROPOSED CONDITIONS 


The heading listing the nodules has output on the terminal. The SCS-TR-S5 
analysis for current conditions is complete. HAP is ready to accept 
changes within the watershed and repeat the SCS-TR-55 analysis. If there 
are no watershed changes to be investigated, proceed to Step Six to exit 
program. 


A. Request the UPDATE module. 

B. As prompted, select the quads requiring changes. 

C. Enter the location and proposed land cover, soil, or slope 
value for each update cell from Tables I, II, and III. In 
this case, only land covers are being changed. 


D. For the given rainfall in STEP FOUR, the updated conditions 
will be routed through TR-55 and the response compared with 
the present conditions. 

E. The user will be prompted for proposed condition maps to 
verify the update input. 

THE TERMINAL INPUTS AND OUTPUTS FOR STEP 5 ARE: 
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STEP SIX 


EXIT HAP 


The heading listing the nodules has output on the teminal. The analyses 
are complete and you are ready to exit from the program by calling the 
•STOP' module. 


THE TERMINAL INPUTS AND OUTPUTS FOR STEP 6 ARE: 
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STOP KKIML IMP EXIT 


These interactive inputs may be placed in a data file prior to HAP 
execution. HAP can be instructed to receive inputs from this file creating 
a quasi-interactive terminal session. This mode is especially advantageous 
when testing the same watershed for many different rainfall rates. 

CONCLUSION 

Because of its interactive, real-time capabilities, the Hydrologic 
Analysis Program (HAP) is very valuable as a tool to assist in decision 
making processes. Although it is optimized for water resource analyses, 
HAP's ability to access any area or cell within the MSDAMP data base and 
then produce statistical tabulations or symbolic maps should make it an 
extremely attractive tool for an array of problems encountered by county 
officials. HAP is currently operational on the MNCPPC's HP3000 mini- 
computer and the U(^CE UNIVAC 1180 main frame computer. 
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